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SUMMARY 


The apanwise lift distributions of vings of low aspect ratio but of 
arbitrary plan form and angle -of -attack distribution have been analyzed 
by two well-established concepts: the virtual -mass concept and the 

Weissinger method. Both concepts have been found to yield the same 
simple integral expression for the spanwise lift distribution in terms 
of the spanwise angle -of -attack distribution. Tables and figures of 
lift distributions of low-aspect -ratio wings have been presented. 

Within certain limitations these distributions are independent of the 
plan form. For plan forms with leading edges swept back behind the Mach 
cone, the results of the analysis are applicable at supersonic speeds. 


INTRODUCTION 


For wings of very low aspect ratio, of the order of 1 or less, at 
small angles of attack the flow may be assumed to be esentlally two- 
dimensional in planes normal to the stream direction. (See reference 1.) 
In the present paper this concept is shown to lead to a simple formula- 
tion of the problem of calculating the spanwise lift distribution of 
wings of very low aspect ratio for various angle -of -attack distributions 
in terms of conditions at the widest section of the wing. The same 
problem is also formulated as a limiting case of the Weissinger 
L-method. (See reference 2.) Some comparisons are drawn between the 
assumptions and limitations of the two approaches . 

Spanwise lift distributions have been calculated for angle-of- 
atfcack distributions vhich correspond to several symmetrical and 
antlsymmetrical twists, as well as to the deflection of flaps and 
ailerons of various spans. The results are presented in a table and in 
several figures. 
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SM30IS 


aspect ratio 

vlng span (at widest section) 
aileron span (total for botb wings) 
flap span (total for both wings) 
chord 




dimensionless chord 

average chord 

local lift coefficient 
wing lift coefficient 

( ^1 / 2 ^ 

lift coefficient on either half of wing 
induced drag coefficient 


rolling-moment coefficient mome nt 

V q.Sb 

bejading -moment coefficient 


lift 

lift on either half of wing 


vertical momentum of virtual mass per unit length 
dynamic pressure 
wing area 
airspeed 

downwash velocity 
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spanwise distance 

dimensionless spanwise 

dimensionless spanwise 
of semi span) 

dimensionless spanwise 
angle of attack 

angle of attack (slope 

angle of attack at tip. 


distance (fraction of semi span) 
distance of center of pressure (fraction 

distance at point of discontinuity in 

of mean-camber surface), radians 
radians 


flap or aileron effectiveness paramister 


( 


5c^/^S \ 

^ci/haj 


circulation 

dimensionless circulation 

/c?c 

loading coefficient 

flap or aileron deflection, radians 
variable of integration corresponding to y* 
trigonometric variable corresponding to y*(cos"l y*) 
variable of integration corresponding to 0 
sweep angle 


or o*oi) 
or |r*) 


ANALYSIS 

The Yirtual-Mass Approach 


According to the theory of reference 1 the disturbance set up by 
wings of low aspect ratio at small angles of attack is essentially two- 
dimensional in planes perpendicular to the free stream. The lift on a 
spanwise strip of width dx may then be obtained from the time rate of 
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change of the vertical momentum of the virtual mass associated with the 
given two-dimensional flow 


dL = M djc 
dt 


= djc 
dx 


= V dM 

•vdiere M is the vertical moinentum of the virtual mass per unit length. 
The lift on the entire wing is then 

^ = If (1) 

Tdiere is the downward vertical momentum at the maxim:jm width of 

the airfoil. Any portion of the wing behind the section of maximum 
width is assumed to be alined with the wake so that it does not affect 
the lift. 

Equation (l) indicates that the lift depends entirely on the 
geometry of the widest section. A physical interpretation of this 
equation is that the wake leaves the airfoil parallel to the surface at 
the section of maximum width and that the vorticity distribution in the 
wake is such as to produce the required downwash distribution in any 
plane perpendicular to the wake. The determination of the spanwise lift 
distribution on the wing then reduces to a problem almost identical in 
form with that of two-dimensional thin-airfoil theory and very similar 
to but simpler than the classical lifting-line theory. 

If the bound circulation on the wing at any spanwise location is 
given by the parameter 



and if the spanwise distance is given by the parameter 


y* = 


y 
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and an associated variatle of integration tj*, then the dowamsh angle 
at any spanwise location y* in the -wake is given by 

--.Jut- (2) 

V 2tcaJ_^ dTj* y* - Ti* ' ' 

From the foregoing considerations this angle must he equal to a, the 
slope of the mean-camber surface at the widest section. If, furthermore, 
the trigonometrical variables 9 and t 3 are defined as 

0 = cos“^ y* 

■B = cos"l q* 


equation (2) becomes 


2itA Jq di3 cos d - cos 0 

Equation (3) may be solved by means of an inversion integral; 


7 = ^ / a sin ■3 log 
Jo 


sin 


3 + 0 


sin 


13 - 01 


d3 




For a constant angle of attack, for instance, equation (4) yields an 
elliptic lift distribution with a lift coefficient equal to ~ Aa in 
agreement with the result of reference 1. 


Three -Quarter-Chord Approach 

Another approach to the problem of calculating spanvise lift dis- 
tributions on wings of low aspect ratio is to consider the limit, as the 
aspect ratio goes to zero, of the equations of the Weissinger L-method. 
This method consists of concentrating the bound vorticity at the quarter- 
chord line, calculating the downwash angle at the three -quarter-chord 



6 


NAGA TW 1973 


line, and setting that do-wnvash angle eq.ua! to the slope of the mean- 
camher surface at the three-q,uarter-chord line • This eq.uation may he 
■written as 

^ ^ ^ ^ . (5) 

•vdiere the dir^nsionless circulation is defined as 

r* = — 
hV 


and the plan -form function F is defined, for tj* ^ 0, as 



and ■vdiere a now represents the slope of the mean-camher surface at 
the three-q.uarter-chord line. 

As the aspect ratio approaches zero, the F function approaches 
the value tan A, ■vdiich is constant if the semispan q,uarter-chord line 
is strai^t. For instance, all plan forms of figure 1 except that of 
figure l(c) have strai^t q.uarter-chord lines. If F is constant the 
second integral in eq^uatlon ( 5 ) vanishes. Furthermore, the dimen- 
sionless chord approaches infinity and its reciprocal vanishes. 

Consequently, the second term of equation ( 5 ) may he expected to 
approach zero quite rapidly as the aspect ratio approaches zero. 



HACA TN 1973 


7 


If the parameter 7, T^ich Is related to the dimensionless 
circulation F* hy 

7 = 

2 

and the trigonometric variables are substituted in equation ( 5 ) and if 
the second term is neglected, equation ( 5 ) becomes identical vith 
equation ( 3 )^ except that a now stands for the angle of attack at the 
three-quarter-chord line. The solution for the circulation is then 
given by equation (i<-) as before. 


EESDIffS 


In order to facilitate the application of equation (4) to cases of 
practical interest, the Integral of that equation has been evaluated 
for several angle -of -attack .distributions -which can be expressed 
analytically as follows: 

(a) The constant angle -of -attack or additional-loading case 
analyzed in reference 1 

The linear antisymmetrlcal-t-virist or damping-in-roll case 
previously analyzed in reference 3 

(c) The symmetrical- and antisymmetrical-t-wist-distrlbution 

cases a = aj.jyx-|°- in -vdilch is the angle of attack at the tip and 

is considered to have equal and opposite values at the t-wo tips in the 
antisyrametrical case 

(d) The flap case in -vdiich a = 058 from j* = to y* = Jq* 

and a = 0 else\diere, agS being the effective angle of attack due to 
flap deflection (some considerations involved in the selection of the 
value of 05 are outlined in the section entitled 'iDlsciission") 

(e) The aileron case in -sdiich a = ±a 5 & outboard of y* = ±yo* 
and a = 0 else-vdiere 

The resulting lift distributions for cases (a) to (e) are presented 
in table 1 and in figures 2 to 5“ All results are presented for unit 
aspect ratio A and unit angle of attack a, -vdiere a is the angle of 
attack at the tip for the t-wist cases and is the effective angle of 
attack due to control deflection a 55 for the flap and aileron cases. 
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The tmst distributions of case (c) may be useful in the analysis 
of structural tvists. If a structural twist distribution is of the 
type considered in case (c), the corresponding lift distribution may be 
taken from the table or the figures directly. Otherwise, a structural 
twist distribution can usually be approximated by the linear super- 
position of two or more of the twist distributions covered in case (c) . 
The lift distribution is then given by a superposition of the lift dis- 
tributions corresponding to those particular twist distributions. 

The results for case (d) may be applied to the case of symtnetrically 
deflected ailerons by subtracting a flap distribution from the constant - 
angle -of -attack distribution for the same angle. Si m ilarly, the results 
for case (e) may be used to calculate the lift distribution due to 
inboard ailerons by subtracting the aileron distribution for the given 
value of yo* from the aileron distribution for 100-percent-semispan 
aileron. 

From the lift distributions the lift, rolling, and bending moments, 
as well as the lateral center of pressure of the wing, may be calculated 
in the form of dimensionless coefficients: 


Cl = J / 7 sin 0 dfi 

2 J 0 


CLi/2 


,«/2 


7 sin 0 d0 


0 


Cj = ^ / 7 sin 20 d0 


0 


1 

Cbm = 2 j 7 sin 20 d0 


y* = 


^^ 1/2 
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The values of these coefficients for the various angle -of -attack 
distributions are given in table 1; the lift coefficient due to flap 
deflection and the rolling moment due to aileron deflection are also 
given in figure 6 as functions of the flap or aileron span ratio. 

The induced drag coefficient may be calculated from the conventional 
lifting- line eq.uation since the drag is not affected by the longitudinal 
distribution of the lift. Consequently, 



a7 sin 0 d0 


(6) 


For the constant angle -of -attack case equation (6) becomes 


Odi = 



as also given in reference 1. 


DISCUSSION 


The tvo approaches to the problem of calculating the span-wise lift 
distributions of vings of low aspect ratio presented in -this paper have 
in common all the limitations of linear potential theory and are subject 
to certain additional inherent limitations as -well. These additional 
limitations and the effect they have on the validity of the results of 
■the t-wo approaches are discussed in the succeeding paragraphs. 

The validity of -the virtual-mass approach has been demonstrated 
for the constant -angle -of -attack case in references 1 and 4 and for 
the linear antis 3 rmmetric angle -of -attack case in reference 5 primarily 
by comparison with more rigorous analytical methods. The validity of 
the Weissinger L-method for plan forms of low aspect ratios has 
similarly been demonstrated for the constant -angle -of -attack case in 
references 2 and k and for the linear antisymmetric case in reference 5 
by comparison -with both analytical and esperimental results. The 
extension of the approach of reference 1 and -bhe specialization of the 
method of reference 2 presented in this paper may therefore also be 
expected to be valid provided certain conditions of aspect ratio, plan 
form, angle of attack, and Mach number are met. 

The aspect ratio has been assimaed to be vanishingly small in the 
analysis. For the purpose of comparison the lift-curve slopes and 
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additional lift distributions calculated by the method of reference 2 
are shown in figure 7 for various plan forms of aspect ratio 1*5 
together with the values calculated hy the low-aspect-ratio analysis of 
the present paper, namely an elliptic lift distribution and a 11ft- 

curve slope of ^(1.5) = 2.356. At an aspect ratio of 1.5 the lift dis- 
tributions for the various plan forms are in fair agreement with the 
elliptic distribution and with each other. The lift -curve slopes of 
the various plan forms are almost the same and are about 15 percent 
lower on the average than that calculated by the simplified analysis. 

This type of agreement does not appear to be sufficient for most 
practical purposes. In reference 1, however, the lift -curve slope is 
shown to be overestimated by only about 14- percent at an aspect ratio 
of 1 as compared with a more exact lifting-surface method. Conseq.uently, 
an aspect ratio of 1 may in general be regarded as the upper limit at 
■vdilch the results of the analysis of this paper may be expected to 
apply. 

The virtual -mass concept utilized in reference 1 does not imply any 
plan-form or angle -of -attack restrictions except that the aspect ratio 
must be low and that there cannot be any abrupt longitudinal changes in 
the span or the slope of the sxirface. As pointed out in reference 1, 
however, it is convenient to apply this concept only to pointed plan 
forms such as those shown in figures l(b), l(d), l(e), l(f)> l(h), l(i), 
and l(j). Wings with re-entrant trailing edges cannot be conveniently 
treated by this method. Furthermore, any part of the wing downstream 
of the widest section must be alined with the widest section. Any flaps 
and ailerons (such as those shown in fig. 1(h)) must themselves be of low 
aspect ratio in order for the low-aspect -ratio analysis to apply to the 
wing. 


The Weisslnger L-method, on the other hand, applies to all plan 
forma except those with curved or discontinuous quarter-chord lines or 
discontinuous chord distributions (it does not, for instance, apply to 
an elliptic plan form) . It has the disadvantage, however, of being 
strictly applicable only to uncambered wings. This observation results 
from the following reasoning: The Weisslnger method is based on the 

slope of the mean-camber surface at the three-quarter-chord line at all 
aspect ratios. For wings of very low aspect ratio the virtual-mass 
concept indicates that the widest section determines the lift on the 
wing. Consequently, the application of the Weisslnger method to wings 
of low aspect ratio can give correct results only if the wing is 
uncambered, so that the slope at the three-quarter-chord line is equal 
to that at the widest section. 

Consequently, the virtual -mass concept is (for practical purposes) 
more restricted as to plan form but less restricted as to camber than 
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the three-quarter-chord method. In view of the agreement of the two 
concepts in the cases to vdiich hoth apply, however, sli^t deviations 
from the plan-form conditions of the one or the camber conditions of the 
other may be expected to have little effect on the results. 

If either concept is used to analyze a wing with deflected flaps 
or ailerons, the question arises as to what value to choose for the 
parameter ag. According to the virtual-mass concept a particular flap 
or aileron angle will be Jiist as effective in producing lift as the same 
angle imposed on the entire section; therefore the parameter 05 is 1 
and the flap or aileron chord ratio does not enter into the problem. 

On the other hand, the three-quarter-chord concept is applicable only if 
the flap or aileron occupies 100 percent of the chord, in \diich case ag 
is again 1. In the more conventional case of partial-chord flaps and 
ailerons it might be expected, by analogy to the traillng-edge cqncept, 
that a value of ag between 1 and the two-dimensional value of 
about 1/2 should be chosen rather than the two-dimensional value for 
the given flap or aileron chord ratio; the determination of the exact 
value would require either experimentation or a more refined theory. 

The angle of attack has been assumed small in the analysis for 
several reasons: At large angles of attack the flows in planes perpen- 

dicular to the direction of flight cannot be considered to be inde- 
pendent of each other, the trailing vorticity is not shed substantially 
in the plane of the wing but more nearly parallel to the direction of 
fli^t, and, finally, at large angles of attack the flow separates and 
potential -flow considerations no longer apply. 

If the leading edge and the sides of the wing (in the case of plan 
forms of the type shown in figs. l(a), l(b), l(e), and l(f), for 
instance) are sharp the flow may detach at very low angles of attack. 

The flow separation affects both the total lift, idiich assumes a 
parabolic variation with angle of attack rather than the usual linear 
one, and the lift distribution. (See reference 6 .) The theoretical 
treatment of this phenomenon transcends the scope of this paper; an 
analysis of a simple case is presented in reference 7 • If the leading 
edge and sides are well rounded, the flow separation and, the resulting 
parabolic variation of the lift coefficient and distortion of the lift 
distribution tend to be postponed to hi^er angles of attack. At lower 
angles of attack the analysis of this paper may be applied to the 
problem. 

The virtual -mass approach has been shown to be valid for supersonic 
speeds (reference l), provided the leading edge is swept back suffi- 
ciently far behind the Mach cone, as well as for subsonic speeds. 
Computation of the downwash in the manner utilized in this paper is also 
valid at supersonic speeds because far downstream of the cross element 
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of a horseslio© vortex the dowawash is the same idiether the vortex Is In 
subsonic or supersonic flow. Consequently, the three-quarter-chord 
concept may also be expected to be valid at supersonic speeds provided 
the leading edge is swept back far behind the Mach cone and provided 
the angle -of -attack, camber, and plan-form conditions are such that the 
analysis of this paper applies. 


COWCniDING EMARZB 


The spanwise lift distributions of wings of low aspect ratio but 
of arbitrary plan form and angle -of -attack distribution have been 
analyzed by two well-established concepts: the virtual-mass concept 

and the Weissinger method. Both concepts have been found to yield the 
same simple integral expression for the spanwise lift distribution in 
terms of the spanwise angle -of -attack distribution. Tables and figures 
of lift distributions of low-aspect -ratio wings have been presented. 
Within certain limitations these distributions are independent of the 
plan form. For plan forms with leading edges swept back far behind 
the Mach cone, the results of the analysis are applicable at supersonic 
speeds . 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Ya., August 25, 19^9 
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Figure 3.— Lift distributions of wings of low aspect ratio for antisyinmetrical twists. 
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Figure 5.— Lift distritutions of wings of low aspect ratio for ailerons of various spans. 
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Figure 6.— Lift coefficient due to flap deflection and rolling-moment 
coefficients due to aileron deflections for wings of low aspect 
ratio. 



